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Graphene, a two-dimensional crystal made of carbon atoms, provides a new and unexpected 
bridge between low and high-energy physics. The field has evolved very fast and very good reviews 
are already available in the literature. Graphene constitutes a condensed matter realization of 
lower dimensional quantum field theory models that were proposed to confront important - still 
unresolved - puzzles of the area: Chiral symmetry breaking and quark confinement. The new 
materials named topological insulators, closely related to graphene, are physical realizations of 
topological field theory. This article reviews some of these topics with the aim of bridging the gap and 
making these condensed matter issues accessible to high energy readers. The electronic interactions 
in the monolayer are analyzed with special emphasis on the recent experimental confirmation of some 
theoretical predictions. The issue of spontaneous chiral symmetry breaking in the model materials 
is also reviewed. Finally we give an extensive description of some recent topological aspects of 
graphene that allow to understand the main aspects of topological insulators. 
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I. INTRODUCTION 

The main conceptual advances in modern physics have usually been prompted by almost simultaneous discoveries 
in different branches. In the past century, statistical physics, quantum field theory (QFT) and condensed matter have 
had their main developments in parallel with the best physicists (Feynman, Einstein, Landau, Wigner) contributing 
to them all. Cosmology and astrophysics develop with inputs from particle physics. This fruitful arena has been 
damaged in the last half of the past century by the fast development and specialization of the different fields and by 
the excess of information and production of scientific papers what has originated a profound gap between different 
areas. The experimental realization of graphene, a two-dimensional crystal made of carbon atoms in 2004 provides 
a new and unexpected bridge between condensed matter and high-energy physics which also involves the fields of 
elasticity, statistical mechanics, and cosmology. 

Graphene is a two-dimensional crystal of carbon atoms arranged in a honeycomb lattice: a single layer of graphite 
[l|. Its synthesis [1, amazing properties [3] and potential applications [1, have granted the 2010 Nobel prize in 
Physics to A. Geim and K. Novoselov. One of the most interesting aspects of the graphene physics from a theoretical 
point of view is the deep and fruitful relation that it has with quantum electrodynamics (QED) and other quantum 
field theory ideas [l|, [^-Q . The connection arises from the fact that the low energy excitations of the system can 
be modeled by the massless Dirac equation in two spatial dimensions. From the QFT point graphene has given 
rise to very interesting developments: the so-called axial anomaly 0] has acquired special relevance in relation 
with the recently discussed topological insulators |9|. |lOl| which provide a condensed matter realization of the axion 
electrodynamics [ll, 12]. Charge fractionalization |l3j| has also been explored in the honeycomb lattice with special 
defects 0, [l^ and quantum field theory in curved space [H, [l3] and cosmological models [H, [ll| have been used to 
study the electronic properties of the curved material. A very interesting development is associated to the generation 
of various types of vector fields coming from the elastic properties or from disorder that couple to the electrons in the 
form of gauge fields [23| ■ 

Many good reviews are already available in the literature [20l - |26j , some of them very theoretically oriented p7l - 
[29} so we will focus on some particular aspects that either have not been sufficiently clarified or are the basis of 
subsequent developments. The topics are also chosen to be in the frontier between condensed matter and high energy 
physics. Under a theoretical point of view the synthesis of graphene has opened a new world where ideas from different 
branches of physics can be confronted and tested in the laboratory. On the electronic point of view it can be shown 
that the low energy excitations of the neutral system obey a massless Dirac equation in two dimensions. This special 
behavior originates on the topology of the honeycomb lattice and has profound implications to the transport and 
optical properties. Although the Fermi velocity is approximately a hundredth of the speed of light, the masslessness 
of the quasiparticles brings the physics to the domain of relativistic quantum mechanics where phenomena like the 
Klein paradox or the Zitterbewegung can be explored. None of these questions arise within the quantum field theory 
approach but its full applicability to the condensed matter system is questionable. It is also interesting to see how 
old low dimensional QFT models, particularly the Nambu-Jona Lasinio and the Gross-Neveu i3l| model, found 
graphene as a physical realization. 

Dirac ferniions appeared in condensed matter before the advent of graphene but are ubiquitous in our days. Again 
here graphene constitutes the cleanest and one of the most robust examples. 

Topology is a invaluable mathematical tool in many areas of physics, particularly, in high energy and condensed 
matter physics. In the latter case we can also find apparently disjoint areas where the dynamics of the relevant 
degrees of freedom in the system are governed by laws based on topological constraints. This is the case for instance, 
of critical systems in two dimensions comprising the XY spin model, dislocations and melting of solids, liquid crystals 
and so on [S^]- Also, the concept of Berry phase will be ubiquitous in the physical properties of crystalline solids, 
as we will see. In basis of its well known properties, we will use graphene as a laboratory to understand most of 
the properties derived from the topology associated to Berry phase. Although the field of topological insulators has 
followed its own path fl^ , some of the most fundamental aspects of the the physical phenomena involved can already 
be found in graphene. Being a very simple material and model, understanding the topological aspects of it opens the 
way to follow the more sophisticated developments. 

The organization of the article is summarized in the contents. 



3 




FIG. 1: Left: The honeycomb lattice is made of two interpenetrating triangular lattices. Right: the dispersion relation. 



II. A SUMMARY OF THE GRAPHENE FEATURES 
A. Monolayer graphene. 

The construction of the free action in condensed matter physics proceeds in a very similar way as in QFT. The non 
interacting Haniihonian is determined by the discrete (crystal) and internal symmetries of the system. The "band 
theory" provides the dispersion relation of the material and the electronic properties for a given electron occupancy. In 
the case of having a metallic system, adding interactions can open a gap and give rise to various non-trivial insulating 
states 1331 or keep the system gapless within the universal class of the Landau Fermi liquid, the standard model for 
metals |34l ]. (A very important exception is the behavior of one dimensional systems that give rise to a different 
universality class, the Luttinger liquid [351 [36|). 

The carbon atom has four external 2s^, orbitals able to form molecular bonds. The crystal structure of graphene 
consists of a planar honeycomb lattice of carbon atoms shown at the left hand side of Fig. [T] In the graphene structure 
the in-plane a bonds are formed from 2s, 2px and 2py orbitals hybridized in a sp^ configuration, while the 2pz orbital, 
perpendicular to the layer remains decoupled. The a bonds give rigidity to the structure and the tt bonds give rise 
to the valence and conduction bands. The exotic electronic properties of graphene are based on the construction of 
a model for the tt electrons sitting in the position of the Honeycomb lattice drawn by the a bonds. Alternatively, 
the mechanical properties involve the a bonds with characteristic energies of the order of 7-10 eV. The low energy 
excitations around the Fermi energy will have characteristic energies ranging from a few meV up to 1 eV. 

Most of the crystal lattices discussed in text books are Bravais lattices and in two dimensions they can be generated 
moving an arbitrary lattice point along two defined vector lattices. This happens in the generalized square and 
triangular lattices in two dimensions. It is easy to see that this is not the case of the Hexagonal lattice. This lattice is 
very special: It has the lowest coordination in two dimensions (three) and it has two atoms per unit cell. As it can be 
seen in Fig. [l] the Hexagonal lattice can be generated by moving two neighboring atoms along the two vectors defining 
a triangular sublattice. This is the first distinctive characteristic that will be responsible for the exotic properties of 
the material. 

The dispersion relation of the Honeycomb lattice based on a simple tight binding calculation is known from the 
early works (37l . [ssj . We will not repeat here the derivation which is very clearly written in any graphene review but 
will instead highlight the main properties. The first is that two atoms per unit cell implies a two dimensional wave 
function to describe the electronic properties of the system. The entries of the wave function are attached to the 
probability amplitude for the electron to be in sublattice A or B. 

The nearest-neighbor tight binding approach reduces the problem to the diagonalization of the one-electron Hamil- 
tonian 

<ij> 

where the sum is over pairs of nearest neighbors atoms i,j on the lattice and a^, are the usual creation and 
annihilation operators. The Bloch trial wave function has to be built as a superposition of the atomic orbitals from 
the two atoms forming the primitive cell: 



(2) 
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The eigenfunctions and eigenvalues of the Hamiltonian are obtained from the equation 



(3) 



where uj is a triad of vectors connecting an A atom with its B nearest neighbors {Si in Fig. [T]), and vj the triad of 
their respective opposites, a is the distance between carbon atoms and e is the 2pz energy level, taken as the origin of 
the energy. The tight binding parameter t estimated to be of the order of 3eV in graphene sets the bandwidth (6eV) 
and is a measure of the kinetic energy of the electrons. The eigenfunctions are determined by the the coefficients Ca 
and Cb solutions of equation The eigenvalues of the equation give the energy levels whose dispersion relation is 



\/3 3 
£'(k) = ±t\j 1 + 4: cos^ —^akx + 4 cos —^akx cos -afcj, , (4) 

in which the two signs define two energy bands which are degenerate at the six vertices K of the hexagonal Brillouin 
zone. The dispersion relation is shown in the right hand side of Fig. [TJ 

The neutral system with one electron per lattice site is at half filling. The Fermi surface consists of six Fermi 
points as can be seen in Fig. [T] (only two are independent). This is the most important aspect of the system 
concerning its unusual properties. The existence of a finite Fermi surface (a Fermi line in two dimensions) in metals 
is at the hart of the Landau Fermi liquid standard model. It implies a finite density of states and the screening of 
the Coulomb interaction. Moreover it allows the construction of the Landau kinematics leading to the possibility of 
superconductivity and other collective excitations in the otherwise free electron system [39j . 

A continuum model can be defined for the low energy excitations around any of the Fermi points, say i^i, by 
expanding the dispersion relation around it: k = Ki + (5k what gives from (jS]) the Hamiltonian 



The limit Vmia^Q "H /a defines the continuum Hamiltonian 

H = vpd.k, (6) 

where a are the Pauli matrices and the parameter vp \s the Fermi velocity of the electrons estimated to he vp ^ 
ita ~ c/300. Hence the low energy excitations of the system are massless, charged spinors in two spatial dimensions 
moving at a speed vp- We must notice that the physical spin of the electrons have been neglected in the analysis, the 
spinorial nature of the wave function has its origin in the sublattice degrees of freedom and is called pseudospin in 
the graphene literature. The same expansion around the other Fermi point gives rise to a time reversed Hamiltonian: 
^2 — vp{—axkx + cfyky). The degeneracy associated to the Fermi points (valleys in the semiconductor's language) is 
taken as a flavor. Together with the real spin the total degeneracy of the system is 4. 

The chiral structure of the spectrum described above and the quantum mechanics nature of the condensed matter 
system (as opposite to QFT), allows to test several predictions of the old relativistic quantum mechanics In 
particular electrons in the graphene system will tunnel with transmission probability one through a step barrier hit 
at normal incidence. This phenomenon known as the Klein paradox (4ll . |43 | has been experimentally confirmed 
(43I I43. A similar phenomenon is the so-called Zitterbewegung or fast trembling motion of the electrons in external 
fields |45l. l46j whose experimental signature has been proposed in (47j . A particularly interesting phenomenon is the 
supercritical atomic collapse [H, [i^l in graphene, a consequence of the large value of the "graphene fine structure 
constant" to be discussed in Sect. IIII A2I Another very interesting phenomenon is the possible realization of the 
Swinger mechanism (soj . i. e. the production of charged particle-antiparticle pairs driven by an external constant 
electric field [Hill. 

The former phenomena are not only curious realizations of relativistic quantum mechanics, they have profound 
consequences on the experimental aspects of the material system most of them not observed before in condensed 
matter. In particular the Klein paradox implies that impurities and other most common sources of disorder will not 
scatter the electrons in graphene. This will evade the Anderson localization [HI] , a very important result establishing 
that any amount of disorder in free electron systems in two space dimensions will localize the electrons turning the 
system to an insulator. The Klein tunneling is also partially responsible for the high mobility at room temperature 
and the excellent metallicity of the system. Zitterbewegung was suggested in [45J as an explanation for the observed 
minimal conductivity of the samples [J] , one of the most interesting aspects of graphene whose origin remains uncertain. 
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It can be proven that the Fermi points are topologically stable against deformations of the Hamiltonian preserving 
the product of the discrete symmetries time reversal and spacial inversion [sl] . The topological stability is based on 
the non-trivial momentum space topology of the Fermi points [s^ and their associated non-trivial topological charge. 
Moreover the wave function around the Fermi points have a non-trivial Berry phase that will be discussed in Sect. 
IVCT] 



B. Multilayer graphene. 



Multilayer graphene consists in a superposition of several graphene layers interconnected by hopping terms. De- 
pending on the relative orientation of the layers several stacking are possible. The two most common, staggered 
(ABA) and rhombohedral (ABC) are sketched in Fig. [2]. The low energy limit of some multilayer systems give rise 
also to interesting QFT models defined in the continuum. Bilayer graphene is one of the best studied cases. It is 
relatively easy to obtain and manipulate experimentally and can be better than the monolayer for some potential 
applications 56|]. 

We will denote the two inequivalent atoms in layer n as (a„, 6„). The simplest model introduces interlayer hoppings 
t only between nearest neighbors. The resulting hamiltonian for bilayer graphene in the vicinity of the Ki Fermi point 
is 



n{k) = 





k* 





t 


k 




















k* 


\t 





k 






and the energy bands are given by 



2 



(7) 



(8) 



In the limit E <^ t, one can obtain an effective hamiltonian for the lowest energy bands [57[. To this end, reorder the 
wavefunctions according (oi, 6i, 02, 62) — ^ (12, bi,ai, &2), so that in the new basis the hamiltonian becomes 



n{k) ^ 




where Hij is a 2 x 2 block. The identity 



shows that, for E t, the substitution H22 — E 
diagonalization of the 2x2 effective hamiltonian 

ye// = 



Hii H12 

H2I E122 



DetiJ-L ~ E) 

2 - E)-^H2i - E^ Det{H22 



(9) 



(10) 



E) 



H22 reduces the computation of the lowest energy bands to the 



H'Ji = i/n - H^2H22'H2i = -- 



/c*2 

e 



(11) 



This effective hamiltonian involves only the atoms (02, &i) which are not linked by t. Under a QFT point of view 
the Hamiltonian is very exotic. It has chirality as the Dirac case but the dispersion relation is quadratic. Under the 
condensed matter point of view it is also a very interesting system where the Fermi surface reduces to a point but it 
has a finite density of states at the Fermi level similarly to standard metals. 

The former analysis can be easily generalized to multilayer graphene with rhombohedral stacking. This type of 
staking includes the links (61 — 02, 62 ~ ■ • • j^at-i ~ o,n) and the effective hamiltonian, which involves only the 
unlinked atoms (ai,&Ar), is given by 



1 





k^ 



*N 







(12) 



It is easy to compute that there is a topological charge associated to this type of multilayer systems around the 
Fermi points whose value is ±N what also guarantees their topological stability Note however that there is 

no topological stability for the Bernal stacking. As we will see in Sect. IIII CI four Fermi interactions give rise to 
various spontaneous symmetry broken gapped phases in perturbation theory, something that does not happen in the 
monolayer system (discussed in Sect. IIII B[) . 
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FIG. 2: (Color online) Sketch of the three possible positions of a given layer with respect to the others in a graphene stack. 
Bernal stacking (1, 2,1,2, ■■■ , is described by two inequivalent planes, while orthorhombic stacking, 1, 2, 3, 1, 2, 3 • • • , requires 
the three inequivalent ones. 



III. QUANTUM FIELD THEORY ASPECTS OF GRAPHENE. 
A. Electron-electron interactions in monolayer graphene 

The effects of the electron-electron interactions in graphene has recently become a hot and highly debated topic. In 
the early times of the graphene boom, it was widely assumed that the main features of graphene could be reasonably 
studied using independent particle models. Theoretical work which did not include interactions where highly suc- 
cessful in explaining the anomalous features observed in the Integer Quantum Hall regime, the dependence on carrier 
concentration of the electron mobility, transport in p-n junctions (the Klein paradox) or the transparency at optical 
frequencies. The lack of interaction effects seemed confirmed by careful measurements of the electron compressibility, 
well explained by independent particle models. The advent of new samples with greatly enhanced homogeneity, and 
with carrier mobilities higher by about two orders of magnitude than the graphene flakes previously used, has changed 
completely the perception about the role of interactions. 

The first hint was the observation of the Fractional Quantum Hall Effect in high mobility suspended samples. This 
phenomenon proved very elusive, but, when found, the associated energy gaps turned out to be higher than in other 
2D electron liquids [Hlfsi 113]. 



1. Long range Coulomb interactions. Graphene versus QED. 

The QFT modeling of Coulomb interactions in graphene has been explained in several review articles [11, [2l|, [23, [6l| . 
We will here summarize the main aspects of it. 

As described in the previous section, the standard non-interacting model for the electronic excitations around a 
single Fermi point in graphene (in units 7i = 1) is given by 

n = VF j Sv4,{Y)-i'd,4'{r) , (13) 

where i = 1,2, V5(r) = il:'^{v)^^, and the gamma matrices can be chosen as = '^2,ly = — cri,7° — . <7i are the 
Pauli matrices and vp the Fermi velocity, the only free parameter in the Hamiltonian. 

The treatment of the Coulomb interaction is what makes a difference between the graphene model and the usual 
QED(2-|-1). In the definition of Coulomb interaction in lower dimensions we must clarify weather we leave in flatland, 
i. e. the electromagnetic field is really defined in two spacial dimensions, or the charges are confined to a plane while 
the Coulomb field propagates in three space dimensions. The scalar potential in classical electromagnetism is defined 
by the Laplace equation in any number of dimensions. Equivalently in QFT it is defined by the lagrangian 

L = Jd'-''+^^x T^.^T^'' , F^, = d^A,-duA^, (14) 

which has two space derivatives. In both cases the effective potential in Fourier space is V{k) ^ 1/fc^. To see the 
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spacial dependence of the interaction we must Fourier transform T^(k) and obtain the know result 

(15) 




In two spacial dimensions the Coulomb interaction in real space is logarithmic. In the real world only the charges are 
confined to the graphene plane while the interaction propagates in three dimensions. This is encoded in the usual 
form of the (instantaneous) Coulomb interaction: 

H.„, = f f ,,, J('.mr^mrf(r.) ,,,, 

27rJ J |ri-r2| 

The interaction (|16p is inadmissible in a QFT approach. It is not only non-relativistic (action-at-a-distance) but it 
is also non-local. We can try to formulate the interacting problem in a QFT language by coupling the fermions to a 
gauge field through the minimal coupling prescription 



Lint =9 J d'xdtf{x,t)A^,{x,t) , (17) 
where g = /Attvf is the dimensionless coupling constant, and the electronic current is defined as 

Now we have to face the problem of coupling the two-dimensional current to a three dimensional gauge potential. In 
ref. this difficulty was solved by using the Feynman gauge 

V^^^ = 0. 

The full Hamiltonian looks like that of (non-relativistic) quantum electrodynamics in two spacial dimensions: 

n^hvF J dh^iv)Y'id^ - ^9Af,)■i'{r) . (18) 

Notice that non-relativistic can be used in two different ways. One is non-covariant in the sense that the time and 
space components of the current have different velocities. This lack of covariance does not invalidate most of the 
rules that make life easy in quantum field theory as the Ward identities. A more severe form of non-relativistic is 
setting the rate vp/c = 0. This amounts to considering that there is no time dependence in the gauge propagator. 
The results obtained in this case can invalidate some of the QFT results as we will see. In ref. Q it was shown using 
renormalization group methods that the hamiltonian ()18|) runs to an infrared fixed point where the Fermi velocity 
equals the speed of life recovering the full Lorentz invariance. 

In any case, despite the formal identification of the hamiltonian ([T8| with QED(2+1), the different photon propaga- 
tors produce different results even if we adopt a relativistic approach (i.e. assume a "retarded" Coulomb interaction) 
from the beginning. 



2. Perturbattve renormalization. 



Ultraviolet divergences arise in quantum field theories due to the singular behavior of the fields at very short 
distances in real space - or at very large energies in Fourier space-. The existence of an underlying lattice in most 
condensed matter theories provides a natural cutoff and the issue of divergences can usually be ignored. Infrared 
divergences appear also in massless theories -and hence there are also present in the graphene - and can be treated 
with different techniques. Here again condensed matter resorts to the finite size of the samples to avoid them. 
Renormalization j62l . l63l | is a prescription to get rid of ultraviolet divergences and construct sensible models where 
physical quantities can be accurately computed. The idea is that the ultraviolet divergences can be cancelled by a 
redefinition of the parameters (mass, coupling constant, wave function) of the theory by adding counter terms to the 
lagrangian. The process is usually done order by order in perturbation theory. If done appropriately, at each order 
one finds finite results independent of the cutoff in the computation of physical observables. 

Depending on the number and badness of the infinities appearing in the computation of the Feynman diagrams 
quantum field theories are classified as nonrenormalizable, strictly renormalizable, and superrenormalizable. In most 
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FIG. 3: Primitively divergent Feynman graphs in QED(4). (a) Fermion self-energy, (b) vertex, (c) photon self-energy. 



cases renormalizable models have an infinite number of divergent graphs -as we go to higher loops in the computation- 
but only a finite "type" of divergences given by the so-called primitively divergent graphs. The models describing 
strong and electroweak interactions in three space dimensions are strictly renormalizable what means that the number 
of primitively divergent graphs exactly equals the number of free parameters in the model. In this case at higher loop 
order we get higher powers of the divergent logarithms and it is possible to perform a formal sum of the perturbative 
series for the parameters. 

Lowering the dimension of the space, strictly renormalizable theories as QED become superrenormalizable: only 
a finite set of graphs need overall counterterms. The number of divergent graphs will still be infinite as if there is 
a one-loop graph G diverging, any graph containing G as a subgraph will also diverge. But if we renormalize G by 
adding a counterterm to it, all these higher order graphs will automatically be finite. The usual rule that assigns 
higher log powers to higher order in perturbation theory simply does not work and there is no problem with it. 

In the case of having a gauge symmetry as in QED, the process of renormalization can interfere with gauge 
invariance and more care is needed to complete the job. If we maintain the postulates of QFT, namely unitarity, 
locality and Lorentz invariance, there are Ward identities that relate some of the renormalization functions. If some 
of the postulates are broken as usually happens in the condensed matter applications, some identities might not work. 

Standard QED(3-f 1) is a strictly renormalizable theory. The coupling constant is dimensionless and, in the massless 
case, the model is scale invariant. The three primitively divergent diagrams shown in Fig. |3] are associated to the 
three free parameters in the theory (the coupling constant and the electron and photon wave functions). The theory is 
strictly renormalizable and all divergences at any order in perturbation theory can be cured by a proper redefinition of 
the parameters. QED(2-|-1) is super-renormalizable. The coupling constant has dimensions of energy what improves 
the convergence of the perturbative series. Massless QED (2-1-1) is ultraviolet finite although it has infrared divergences 
[g^ . Graphene sits in between QED(3) and QED(4). Of the three graphs shown in Fig. [3]only the electron self-energy 
diverges and in the case of considering a static photon propagator only the spatial part has a logarithmic divergence 
leading to the renormalization of the velocity parameter in (??). In the static model widely used in condensed matter 
applications, the electron and photon propagators in momentum space are given by 

Go(fc ,k)=^-^^^^^^^-^, (19) 

no(k) = i^. (20) 

The renormalization functions associated to the electron self-energy are defined as: 

Go 1 - S(fcO,k) = Z;i/'(fc",k)[fc070 - Z,(fc",k)z;7 • k], (21) 

The extra parameter v appearing in (jl9[) in the graphene case has an associated renormalization function Zy which is 
a new feature characteristic of graphene. In the Lorentz invariant massless QED in any dimensions, only the electron 
wave function renormalization is associated to the the electron self energy. 
The wave function renormalization 

as(w,k) 

Z-i> ^ |a;=0 , (22) 

defines the anomalous dimension of the field 

^ dlogZ^/dl, (23) 

{I is the RG parameter that for a hard cutoff A is ^ = log A) and, hence the asymptotic behavior of the fermion 
propagator: 

G(w,k)~„^o^. (24) 
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The function Z^, defines the running of the couphng v whose beta function is defined as 

/3„ = dZ^ldl. (25) 

The fixed points of the system are determined by the zeros of this beta function and the running of the velocity is 
given by the local behavior of the beta function around the fixed point ^6^] . 

2 

It can be seen that in the static model the perturbative series is organized in terms of a coupling constant ac = ■ 
Plugging- in the Fermi velocity measured in nanotube experiments vp ^ c/300 [65], and the bare electron charge, we get 
aQ ~ 2.3 — 2.5. Although this value puts the model in the strong coupling regime, the results obtained in perturbation 
theory agree with the experimental findings. A possible explanation lies on the fact that the effective constant in a 
material has an extra parameter: the dielectric constant taken as one in the vacuum. In samples on a substrate the 

2 

dielectric constant reduces the effective coupling constant: olq = 4^egt,p ; where eg includes intrinsic contributions and 
effects due to the environment in which graphene is immersed. A simple estimate gives eg = {£air + £subs)/^ ~ 2 — 3 
for typical substrates. The actual value of the intrinsic contribution to the dielectric constant due to the graphene 
layer itself is still object of vivid controversies. By extrapolating measurements of the excitation spectra in graphite a 
large intrinsic dielectric constant in single layer graphene of the order of sg ~ 13 was proposed in 66]. An independent 
estimate of a in graphene has been obtained from measurements of the carrier-plasmon interaction in samples with a 
finite carrier concentrationig^]- The result, a « 2.2 is consistent with previous analysis \6§\ and from recent numerical 
calculations [g^. Theoretical values of the order of eg ~ 4 have also been reported using various approximations 

MM- 



3. Velocity renormalization. 
The computation of the graph in Fig. [3] (a) gives 

E«(k)^|z;7-klog^, (26) 

where A is a high energy cutoff and g is the bare coupling. The logarithmic divergence can be fixed by redefining the 
fermion velocity v which becomes, after finishing the renormalization procedure, energy-dependent. The onle loop 
renormalization group equation can be solved to obtain 

— 7T7T = 1 - -rr- ^og(E/Eo), (27) 
vp(£jo) 167r 

which relates the value of the Fermi velocity at an energy E, with that at a reference energy Eq, assuming that the two 
energies are sufficiently closed. Eq is not a cutoff of the order of the bandwidth but any energy where the Fermi velocity 
is determined by experiments. The corresponding diagram in massless QED(3) does not have a logarithmic divergence 
because the photon propagator has an extra inverse power of momentum. In QED(4) this diagram induces a wave 
function renormalization. It is interesting to note that the electron charge in the graphene model is not renormalized 
(after fixing the velocity divergences, the photon self-energy is finite at all orders in perturbation theory) but the 
velocity renormalization induces a renormalization of the graphene structure constant similar as the one obtained in 
QED(4). This issue will be explored further when discussing the experimental observations in Section llll A 41 

Since the bare coupling constant of graphene can be large at the experimentally accessible energies, the former 
renormalization scheme can be improved by performing a 1/N expansion [t^ ] where N is the number of fermionic 
species. In the case of graphene with the spin and Fermi points degeneracies the phy sical value is = 4. This 
procedure was followed in the early publications [tI, [tI] and was later retaken in [75l - l77f . 

The simplest non-perturbative calculation amounts to compute the electron self-energy graph in Fig. [3l[a) with an 
effective propagator obtained from the resummation of the planar diagrams dominant in the 1/N approximation. The 
IIoo component of the one loop photon self-energy in Fig. [Sfb) is given by 

hence, the effective Coulomb potential obtained by summing the geometric Dyson series for the photon self-energy 
in the so-called random phase approximation (RPA) is 



(29) 
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The electron self-energy computed with this expression gives the following p function for the Fermi velocity [TJ, |75| : 

where I is the scaling parameter [l = log A in a cutoff scheme) , and g = ^^^^ . The large N limit amounts to take the 
limit TV — oo keeping g fixed. As we see the dependence on g is non-perturbative and the growth of the velocity at 



low energies is slightly different than that in ((27)) . 

A remark is here in order. In graphene many classes of lattice defects can be described by gauge fields coupled 
to the two dimensional Dirac equation l2^. The standard techniques of disordered electrons fw\ can be applied to 
graphene by averaging over the random effective gauge fields. A random distribution of defects leads to a random 
gauge field, with variance related to the type of defect and its concentration. These random fields when treated 
perturbatively, lead to a renormalization of the Fermi velocity which makes it to decrease at lower energies opposing 
the upward renormalization induced by the long range Coulomb interaction. The simultaneous presence of interaction 
and disorder gives rise to new interesting fixed points. The issue was analyzed in [t^. The most interesting case arises 
when considering a random gauge potential which models elastic distortions and some topological defects. There is a 
line of fixed points with Luttinger-like behavior for each disorder correlation strength A given by v*p = 2eVA. An 
extensive analysis of the issue in the large N limit is done in [75]. 

Another interesting prediction done within this framework is the linear dependence of the quasiparticles lifetime 
with the energy [73l | , a distinctive of the marginal Fermi liquid behavior [80] . Experimental evidences of this behavior 
have been described in I8ll-l84l. 



4-. Experimental confirmation of the Fermi velocity renormalization. 

The RG prediction of the velocity renormalization has been very recent confirmed in several experimental reports 
[ssl - fSTj following earlier more indirect evidences [H, H^. The experimental evidence more directly related to the 
physics discussed in this review is that in [s^. The experiment measures the effective mass of graphene at different 
carrier densities in high mobility suspended samples and also in graphene on a BN substrate. The experiment probes 
the suspended graphene samples at low energies in a range form lOOmeV to 0.2 meV never reached before. They use 
the temperature dependence of Shubnikov-de Haas oscillations to infer the dependence of the Fermi energy, e^^, on 
the area of the Fermi surface, Sf- The effective mass is defined as 

dSp 



so that for graphene 



= ^ (32) 



A comparison of experimental results an fits based in eq. (|30l) are shown in the left hand side of Fig. ID 

The behavior of the inverse running coupling constant in graphene taken from the data in [85] is shown in the right 
hand side of Fig. |4l 

The difficulty of achieving these results is similar to that of going several orders of magnitude higher in energy to 
confirm the fine structure constant renormalization [6l[. It is interesting to compare with the case of QED(3+1). 
The fine structure constant is defined by uqed — e^/ (47reoc) where e is the electron charge, c the speed of light, and 
eo the dielectric constant of the vacuum. As such, it is a dimensionless quantity whose universally known value is 
= 137. From its definition, measuring the fine structure constant amounts to measuri ng th e electron charge with 
the higher possible accuracy since the other quantities in the definition are real constants [l8ll |. 

While the value of e obtained through Hall resistance measurements reached an accuracy of six significant digits 
[ooj . the actual value got from the magnetic moment of the electron has an accuracy of 0.7 parts per billion [oij . 
This allows a determination of uqed of similar accuracy. The calculation includes contributions from 891 Feynman 
diagrams and is one of the most demanding comparisons of any calculation and experiment ever performed [92] . The 
impressive accuracy reached at the solid state energies (around lOOmK considered as zero for the high energy running 
coupling constant) is spoiled at energies of the order of the proton mass (IGeV) where strong interaction diagrams 
contribute to the determination of uqed ■ A relatively recent measurement done at the Large ElectronPositron collider 
(LEP) at CERN provides a value of a~i(M|) = 128.9460.030 [H, H^. From this we can see that the slope of the 
variation in two order of magnitude is similar in the two cases. 
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FIG. 4: Left: Experimental measurements of the effective mass of carriers in high mobility graphene samples. Red dots give 
results in suspended samples, and blue dots are for samples on BN substrates. The red and blue lines are fits using eq. |30] 
Red: Eg = 1. Blue: so = 4. The dotted orange line shows a fit obtained neglecting the renormalization of the Fermi velocity. 
Right: Plot of the Fermi velocity renormalization from the measurements done in [8a |. 



5. Retarded Coulomb interactions and the structure of the perturbative series. 

The success of the effective description of the long range Coulomb interaction in graphene by means of a static 
Coulomb potential described before that lead to the observation of the Fermi velocity renormalization means that 
retardation effects of the order of vp/c ^ 1/300 are ineffective in the today's experimental accessible region. Nev- 
ertheless there are important issues concerning the structure of the QFT model of graphene that we will address 
here. 

First the inclusion of a covariant photon propagator 

instead of ([20]) leads to important conceptual differences in the infrared behavior of the model. 
The one loop computation of the electron propagator in this case gives the result 0: 

1 1 - + 1 - 4i;2 

^ Fi 2\V2 arccos v + 5- , 34) 

where the speed of light has been put to one. The equation 13^ = Q has the solution v = 1 which in our units means 
that we have a fixed point of RG where the electronic velocity equals the speed of light. 

Since the electric charge e is not renormalized, the value of the coupling constant at the fixed point g* is the fine 
structure constant of QED: 

g*=aQED^^- (35) 

This issue was initially studied in [8] was later confirmed in a non-perturbative approach in [95| . Moreover there is 
a wave function renormalization at the one loop level that has been may give rise to a non-trivial critical exponent. 



B. Chiral symmetry breaking or gap opening in monolayer graphene. 

Opening a gap (or generating a mass) in graphene is crucial for the electronic applications. Under the QFT point 
of view, the electron mass (gap) is protected by the 3D version of chiral symmetry and hence a gap will not open by 
radiative corrections at any order in perturbation theory. As discussed in |54| interactions or disorder can open a gap 
in the sample provided that the product of time reversal times inversion symmetry is broken (staggered potential, 
external magnetic field), or when the two Fermi points are involved (Kekul'e distortion). A discussion of the various 
" masses" that can be generated in graphene when the valley and spin degrees of freedom are taken into account will 
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be done in Sect. IV El when discussing the topological insulator aspects of graphene. Even if the two Fermi points 
are taken into account in a four dimensional formalism, that a gap will be spontaneously generated by interactions in 
monolayer graphene is very unlikely and there are so far no experimental evidences for it. 



1. Long range Coulomb interactions. 



One of the interesting features of QED(3) is that a fermion mass can be generated dynamically, breaking chiral 
symmetry. Starting with massless bare fermions, they can acquire a dynamical mass through nonperturbative effects. 
Writing the full fermion propagator as 

a nonzero solution for B(p) implies a nonzero condensate, and signals dynamical mass generation. The infrared value 
of the dynamical mass function defined by 

m(0) = B{Q)/A{Q) can also be used as order parameter. 

The standard non-perturbative approach to study dynamical fermion mass generation is to solve the Dyson- 
Schwinger (DS) equation: 

^j^T,ip,k)D^,ip-k), (36) 

where T^{p,k) is the vertex function and Dfj^i, the photon propagator. Eq. p6|) can be decomposed in a couple of 
equations for A(fc^) and B{k^): 

1 f (fik 

^^P^^ = 1 + V y j^tr[j.kj,Gik)T,D^,{p - k)]. (37) 



1 r rPk 

B{p') = —J j^tr[^,G{k)V.D,,{p-k)]. (38) 

If the DS equation for B{p2) has only vanishing solution, the fermions remain massless and are stable against gauge 
fluctuations. 

The issue of chiral symmetry breaking (CSB) in massless fermion models in QFT is far from being settled. Lower 
dimcnsion a.1 QE D, in particular QED(3) was proposed in the early times as a toy model for CSB 96 99] and confine- 
ment [lOOl . Il0l| . It is clear that no mass will be generated in the weak coupling limit. But, as we described in Sect. 
nil A 21 the bare coupling constant of graphene is not small. 

In the strong coupling regime the best studied approximation is based on a large number of fermionic species, 
the 1/Nf expansion (not to be confused with the 1/N methods used in QCD where N is the number of colors of 
an SU{N) gauge theory). The crucial difference between the perturbative and non-perturbative approaches is that 
the DS equation is nonlinear which makes it possible for chiral phase transition to happen at the bifurcation point. 
Spontaneous chiral symmetry breaking in QED (3) in the 1/N approximation usually requires an unphysical number 
of fermionic species N < 2 — 3.3 102]. The situation with graphene is even worse since there is no mass parameter 
to begin with. The gauge propagator is crucial in these approaches and hence it is important to distinguish between 
QED(3) and graphene. This problem has been addressed in jl03j with the conclusion that for small g there is no 
solution wit h sp ontaneous chiral symmetry breaking. A very good updated review with a fair list of references can 
be found in [lOj]. 

The simplest non-perturbative calculation that can be done to study the issue of the gap opening in graphene is the 
RPA type described in Sect. IIII A3I [74]. The absence of a constant term in the inverse electron propagator ensures 
that no ma ss is generated in this approximation. The 1/N expansion approach to the problem has been revised 
recently in [l05| with the conclusion that a gap will not open for the physical values of the electronic degeneracy in 
graphene (N=4) . A variational approach to the excitonic phase transition in graphene including the renormalization 
of the Fer mi ve locity [l06j | also produces quite negative results. At finite temperature the critical Nc is temperature 
dependent jl07j while a finite chemical potential leads to strong suppression of the critical fermion flavor Nc and the 
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dynamical fermion mass in the symmetry broken phase |108j . Volume effects have been analyzed in Il09j trying to 
explain the discrepancies between the continuum and lattice results. 

The situation can be improved by the presence of an external mag netic field given rise to the so-called the magnetic 
catalysis [ll3|- This possibility was put forward in the early papers Illl4ll3 | and is discussed at length in the review 



2. Four fermi local interactions. 

Short range interactions, such an onsite Hubbard term U are irrelevant in an RG sense [61} in the weak limit. This 
is due to the vanishing density of states at the Fermi level that occurs at Fermi points, a characteristic of graphene. 
This is a very important aspect since short range interactions are responsible of the Fermi liquid properties of usual 
metals (having a finite Fermi surface) in two space dimensions. Nevertheless the density of states at low energies can 
be increased by the presence of disorder and a finite temperature what in turn, enhances the effect of short range 
interactions |ll 4|. These interactions can be relevant in the strong coupling regime of the hexagonal electronic and 



optical lattices 115l [TT 



Four Fermi interactions were considered in low dimensional models in high energy physics again searching for chiral 
symmetry breaking and confinement. One of the most popular models is the Nambu-Jona Lasinio (NJL) in two space 
dimensions [30] . The original Nambu- Jona-Lasinio model was an extension of the BCS-theory of superconductivity 
to the domain of spontaneous symmetry breaking in the strong interaction, one of the best examples of hybridization 
between condensed matter and high energy physics. The lagrangian density is given by 

L = -i^j^'d^ + g [(§'5)2 - (§75*)2] , (39) 

The four-fermion interaction is attractive for fermions and antifermions of opposite chirality. Fermion-antifermion 
pairs form bound states, which form a condensate: 

< >7^ 

that changes the nature of the vacuum of the theory. 

In the condensed matter language the NJL model is similar to the Hubbard model that has been extensively studied 
in lattice models. T his type of interactions has been used in the Honeycomb lattice searching for the possibility of 
gap opening in [t^ Ill7l4l21 [. More recently lattice gauge field theory techniques are being used to explore this 



possibility. An interesting recent analysis of the possible phases arisen in the Jona-Lasinio model in three dimensions 
is do ne in [122| |. Non-perturbative constructive field theories techniques have also been used and are described in 
[l23| . Although the issue remains controversial, it is fair to say that there is no conclusive evidence up to now that 
the physical parameters of graphene lie in the region where chiral symmetry breaking will occur. 

A gap can open in the sample by various extrinsic means although without altering the basic properties of graphene, 
the produced gaps are very small and not entirely controllable. Small gaps hav e be en experimentally reported in 



samples on a substrate whose lattice is commensurate with that of graphene 124| . Graphane, a hydrogenated 
compou nd ha ving hydrogen atoms forming sp^ bonds has a reasonable gap and is an interesting material on its own 
grounds [l25j . 



C. Electronic interactions in multilayer systems. 



The situation concerning electron-electron interactions is much more interesting in bilayer graphene. The quadratic 
dispersion relation and a finite density os states at the Fermi surface enhances the role of four Fermi iriteraction and 
a variety of broken symmetry phases can arise similarly to what happens in the one dimensional case [126] . 

New experiments in high mobility samples hint to the existence of novel phases at zero carrier concentration and 
low temperatures. The first indication came from the electron compressibility measurements reported in |l27l |. These 
experiments showed a rise in resistivity near the neutrality point, consistent with a tendency towards an insulating 
state, although the resistivity never went above a value of a few thousand ohms. The existence of a gapped insulating 
phase was reinforced by extrapolations from measurements made at finite magnetic fields[128], which suggest a gap 
of about A « 2 meV. 

The experiments reported in [l29j | present a different picture. The density of states at the Fermi energy is inferred 
from careful measurements of the carrier density and temperature dependence of the resistivity in a number of 
high mobility suspended samples. The results show a crossover to a low temperature regime where the density 
of states is significantly reduced. The band structure of bilayer graphene changes from parabolic to four Dirac 
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cones, due to trigonal warping. The crossover found experimentally occurs at an energy of about 6 meV, which is 
substantially larger than the crossover energy related to trigonal warping, aboul 1 meV. Hence, the results suggest a 
spontaneous symmetry breaking associated to interactions. This explanation is supported by measurements performed 
at finite magnetic fields, which imply that the lowest Landau level is fourfold degenerate, while independent electron 
calculations give an eightfold degenerate n = Landau level. In any case, the results reported in (129j do not show a 
finite gap at any concentration. 

A last batch of recent experiments is discussed in[T30l|. The conductivity of s suspended sample at the neutrality 
point was measured as function of bias voltage, magnetic field, and perpendicular electric field. The results indicate 
a gap of for bias voltages between abut -3 and 3 meV at zero magnetic field, zero electric field and zero carrier 
concentration. The differential conductance shows peaks at voltages above and below this gap, while they tend 
towards a constant value at larger bias voltages. Peaks adjacent to a gap are a hallmark of tunneling experiments into 
superconductors, although there is no apparent reason for these peaks to show up in a d. c. transport experiment. 

All the experiments described above coincide in that a broken symmetry phase due to electron electron interactions 
seems likely in bilayer graphene near the neutrality point and in the absence of a perpendicular elec tric field. They 
differ in many details, however, and the experiments in[l29j suggest a gapless phase, while those in jl27j |. and especially 
those in[13fl] seem to imply the existence of a gap. Th e exp eriments anal yze suspen ded samples with hig h ele ctron 
mobility, although the mobilities differ, being highest in |l29| . and lowest in |l27l . Il2g |, while t he samples in [l30f show 
an intermediate value. The experimental setups also differ: the measurement in 127| , Il28j give the values of the 
compressibility in samples with one back gate, while the experiments in 129] and in 130l | report the d.c. conductivity 
in samples with no and with two gates respectively. 

On the theory side, it was soon realized that the combination of parabolic bands and short range, screened in- 
teractions lead to logarithmically divergent susceptibilities. The variety of electronic degrees of freedom in bilayer 
graphene allows fo r ma.ny p ossible broken symmetry ground states, which, in turn, can be induce d by appropriately 
tuned interactions jl3ll - |l39l ]. There is a ga pless , nematic phase consistent with the experimen ts iiijl29l |. and different 
gapped phases could explain the results in |l27l |. while the interpretation of the experiments in [l30j | remains less clear. 
Some of the proposed phases break time reversal symmetry, leading to states with similar properties to a 2D elec- 
tron gas in the Integer Quantum Hall regime 0- Recent c alcu lations suggest that many possible phases are almost 
degenerate, with energy differences per atom below 1 meV |l4Cll |. 

Interactions in a graphene bilayer are screened, do that they decay faster than the inverse of the separation beyond 
some distance. The nature of the broken symmetry phase may depend on the value of the screening length[141J. This 
leads to the intriguing possibility that the two phases which look more stable, a gapped layer antiferromagnet and 
a metallic nematic phase, could have been observed in different experiments, depending on the number and position 
of gates. Further complications are introduced by stra ins, which are probably unavoidable in suspended samples. 



Disorder can induce, in some circumstances, local gaps[142[. This situation is reminiscent of other materials with 



many competing interactions, such as the cuprate superconductors or the manganites, where the interpretation of the 
low temperature phases is still debated. At least, graphene has a simple, stochiometric composition. 



IV. GAUGE FIELDS FROM LATTICE DEFORMATIONS. 



One of the most fascinating facts about graphene is the possibility to find experimental realizations quite abstract 
theoretical ideas. There is a particular aspect that resembles what happened in the physics of the standard model 
in the last century: Gauge bosons were postulated on theoretical necessities and found experimentally afterwards [] . 
The emergence of gauge fields associated to lattice deformation has been discussed at length in ref. [20]. We will 
here summarize the main aspects before revising the experimental observation of them that has taken place in the 
last year. As described in [20| gauge fields have emerged from three different sources: Topological defects, elastic 
distortions, and from a covariant, geometrical formulation. The physical origin of them is rooted in the tight relation 
between the spinorial nature of the carriers in graphene and the lattice structure. As described in the introduction 
the spinor degree of freedom is associated with the non trivial geometry of the Honeycomb lattice. 

ARPES measurements show that, despite significant deviation from planarity of the crystal, the electronic structure 
of exfoliated suspended graphene is nearly that of ideal, undoped graphene |84| . 



V. TOPOLOGICAL ASPECTS. 



Here we will review some aspects of the topology the momentum space in fermionic systems on a lattice. We will 
use graphene and the Quantum Hall effect as test systems to see these ideas in action. 
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A. Topology and Condensed Matter Physics. 



In a very broad sense, topology is the branch of mathematics devoted to the quahtative study of forms and their 
classification in terms of the invariance under certain transformations. For instance we can globally classify three 
dimensional surfaces by counting the number of holes they have. It is very well known that the surface of a sphere 
is topologically equivalent to the surface of a bottle whereas a donut is topologically equivalent to a mug or a pipe. 
In condensed matter physics this classification of objects in terms of shapes is useful although not so obvious as the 
examples mentioned before. During the last century, the classification of the quantum phases of matter has been 
done based on three major pillars: the band theory of solids, the Landau theory of the Fermi liquid, and the concept 
of spontaneous symmetry breaking. Topology plays a role in the third pillar in terms of the notion of the order 
parameter and homotopy theory. It is not our intention to develop these ideas here and we refer to the literature |l44l | 
for a general and complete overview of the problem. We will nevertheless mention some of the most the prominent 
examples: The planar vector model whose isolated vortices in the vector order parameter cannot be continuously 
deformed to any smooth spin configuration in the plane. The liquid helium where the vortices appear in the phase of 
the superconducting wavefunction. The theory of defects in crystals |l45j | and so on. 



B. Quantum Hall effect in graphene. 

As it was discussed in Sect. |lT]the bandstructure of graphene possess two distinctive features when compared with 
other realizations of a two dimensional electron gas: the linear shape at low energies and the presence of a SU (2) 
quantum number different from the spin, often referred in the literature as the valley spin or pseudo-spin. These two 
distinctions of graphene are the key features when we consider the topological properties of this system related to the 
Berry phase. 

Let us consider a single layer of graphene under the effect of an external homogeneous magnetic field, described by 
the vector potential A — B{—y, 0,0). In the basic description of the Quantum Hall Effect (QHE) in graphene we do 
not need to use the full description in terms of the two species (valleys) of Dirac fermions just because the external 
magnetic field breaks time reversal symmetry and both species will contribute equally to the Hall conductivity. We 
will extensively use of the so called magnetic length I b defined as Ib — ^^7= • 

The low energy hamiltonian for one valley in graphene can be written as (h ~ 1) 

H ^ vcrU = V {a^n^ + ayUy) , (40) 

with Hj — —idj + eAj. In terms of the ladder opperators a = J^(n^ — iHy) and — j^i^x + ^n^), satisfiying 
[a,a^] = 1, the hamiltonian H reads 

It is important to note that here we are considering the situation of a perfect crystalline layer of graphene with 
inversion symmetry so there is no term proportional to u^. The presence of such important term and similar ones 
will be considered later on. The Schrodinger equation with the hamiltonian (j4ip can be easily solved in terms of the 
solutions of the harmonic oscillator. The important difference here with the standard two dimensional electron gas is 
that the hamiltonian (j41[) is a two dimensional matrix, and the eigenvectors are two dimensional spinors of the form 
V'n = (""nj Vn)'^ ■ In terms of the second component of this spinor, Vn, and the number operator fi = a^a, we conclude 
that ^S-nVn — E^Vn — nvn, with n being an integer. Without entering in more details, we see that for each value of 

n we get two eigenvalues E±_n = ±u-y/2|n|/ZB = ±Wc-\/|n[. An important observation from this equation is that there 
is an eige nstate corresponding to the value n = with no analogue in the standard non relativistic two dimensional 
case pq . 

In order to extract the information concernin g th e topology of this system we shall calculate the effective action 
through the use of the Landau levels structure [147| . We will consider the situation of zero temperature and finite 
chemical potential fi. Within a path-integral approach, the effective action is calculated through Seff = —ilnZ. 
After integrating out the fermion fields in Z the effective action reads 

Seff= / ^ E ^^{^-^^-En). (42) 

n— — 00 
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Because En does not depend on the momentum k each Landau level is (highly) degenerate. The degeneracy of each 
Landau level was first stimated by Landau himself to be -^f^^- By adding an small parameter 77 — >■ 0~^sign{uj) in the 
argument of the logarithm, the integration over the uj variable can be performed, leading to the following effective 
lagrangean 



leBl 



n— — 00 

or, splitting all the contributions from n = and n ^ we get 



^eff = ^Ia^I + 1^ E (l^ + ^-^1 + 1'^ - ^cVH\) ■ (44) 



The averaged particle number can be calculated from (j44p by taking the derivative of Cef / with respect of its conjugate 
variable /i: 

^1^- ^Er=i(^*WM + -c/H) (45) 



Before entering in the details of cq. pS]) let us see what is the actual meaning of the magnitude rt as a function 
of the magnetic field B. From (|46p we see that all the terms are proportional to B (without loss of generality we 
shall assume that B > 0): n — anB, being an the Hall conductivity. We know that the average particle number 
is nothing but the temporal component of the gauge invariant electronic current; n = , and that the magnetic 
field _B, being the z— component of the vector magnetic field B, can be written in terms of the vector potential v4^, 
B — — d2Ai — diA2. Using these considerations and noticing that the system is gauge invariant, the previous 
relation n = uhB can be written in a the following gauge invariant way: 

J'' = cjHe^^'Pd.Ap. (46) 

The equation of motion PS|) which describes the effective dynamics of the electrons under the effect of an external 
magnetic field can be derived from the well known Chern-Simons action: 

^^cs = ^ y d^xe^-'PA^d^Ap - Jf^Ap. (47) 

So from the knowledge of the Landau levels, we have arrived to the well established effective lagrangean describing 
the physics of the integer QHE. The behaviour of the Hall conductivity an can be obtained by inspection of eq. pS)) 
and it is plotted in fig.©. The first two terms in the right hand side of ((46)) have the form of a staircase changing its 
value by integers when the chemical potential crosses the energy of each Landau level. This is exactly the behaviour 
found in the standard two dimensional electron gas, plotted also in fig®. The difference is thus the last term, 
sz.gn(/i), wich adds a remarkable change of one half to the entire Hall conductivity. Being remarkable, this is actually 
not in odds with the result of Thouless and others, who showed that an must take integer values. The solution of 
the paradox is that we actually have two Dirac points in the bandstructure of graphene (garanteed by the Nielsen- 
Ninomiya theorem) and we have to multiply the entire result by two Dirac species and by two spin orientations, so 
the actual values taken by the Hall conductivity are [N e Z) 

CJH = 2-{2N+1). (48) 



C. The Callan-Harvey mechanism and the edge states in the QHE. 



We will comment on two aspects concerning the effective Chern Simons lagrangean. One quite general, related to 
the structure of the effective CS lagrangean (I47p and the existence of chiral states at the boundaries of the system. 
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FIG. 5: Hall conductivities for graphene in the single cone approximation (continuous red line) and the two dimensional electron 
gas (blue dashed line) as a function of the filling factor v. The (green) dotted dashed line corresponds to the case of graphene 
with a mass. 



and the second, more specific to graphene, regarding the appearance of the parity anomaly when a mass term is added 
to the hamiltonian (l40l) . 

It is easy to check that when dealing with ideally infinite systems, the laction defined in (|47|) is gauge invariant 
under the gauge transformation — >■ + 9^A, as long as A is some well behaved function in the entire space. 
However, when a bounded region is considered (or more generally, a domain wall), the action (|47p is no longer gauge 
invariant: 



S'cs = Scs-^ f d^xA{x)e''''Ph^{x)d,Ap, (49) 
47r Jon 

where fif^{x) stands for a unit vector perpendicular to the boundary at any point. Without loss of generality, we can 
choose the boundary to lie on the line y ~ and thus the vector to point in the (0, 0, 1) direction. Then, the 
second term in the right hand side of P51) becomes (i, j = 0, 1) 



SS^-^ cPxK{xy^d,Aj = / £xK(x)€'^F,j. (50) 
47r J in J 

For a generic regular function A inside the bounded region fi, the right hand side of (I50p is nothing but the expression 
for the (1 + 1) chiral anomaly, so through the noninvariance of the Chern Simons action in a region with boundaries, 
we reach the interesting conclusion that we can recover the gauge invariance for the system in a bounded region if we 
add to the system a chiral fermion whose anomalous contribution to the action exactly cancels the noninvariant piece 
of ((50|) . So through this mechanism, known as the Callan-Harvey mechanism |148|, we can rationalize the existence of 
chiral edge states at the boundaries of a system in the quantum Hall regime. Because an is an integer, we see that 
the value of the quanticed Hall conductance tell us the number of chiral edge states we have in our system. Let us 
stress that this argument lies only on the gauge noninvariance of the Chern Simons action and is applicable to any 
two dimensional fcrmionic system, graphene or the standard two dimensional electron gas. 



1. Parity anomaly and QHE. 

In Quantum Electrodynamics in 2 + 1 dimensions an interesting phenomenon emerges when a single species of 
massive Dirac fcrmions in the continuum is considered pj^]. The action for such a system coupled to an abelian gauge 
field A^ is (again, in units h = c = 1) 



S = / d-^xV^ (7^(z9^ + eA^) - m) ^. (51) 
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It is widely well known that the mass term in the action (1511) breaks time reversal symmetry and parity just because 
in order to be symmetric under any of these two symmetries, m must change to — m, fact that only happens for m — 0. 
This lack of invariance under these two discrete symmetries will have as a consequence the emergence of a Chern 
Simons term in the one loop effective action for the field A^. 

The one loop contribution to the polarization function is given by the bubble diagram, and it takes the form, in 
terms of the Feynman propagator for the Dirac field: 



U^-'ik) = e^Tr [ f^r—, Y ^ • (52) 

J Stt-^ jP[qp + kp)~m ■y'^q„-m 

Instead of evaluating ([5^ in its full glory, let us take a short cut to see how the Chern Simons term emerges. We 
see from (|47p that the Chern Simons term is linear in derivatives, so in terms of momenta, it will be linear in k. So 
We can expand ([52l) in powes of k and calculate the linear term. Differenciating ([52]) and setting A; = it is found 
that the coefhcient is given by the following integral 



e'Tr / r Y ■ (53) 

After performing the trace, among all the terms in ([55)) there is one proportional to m^e'^'^P: 



This integral can be easily evaluated taking into account that it depends on and it is insensitive to the sign of the 
mass (also, = m?m = |mpm). The final result is that the effective action for the gauge field contains a Chern 
Simons term through the one loop contribution to the propagator: 



m 

An \m\ 



(55) 



We must stress that this expression has been calculated after expanding the polarization bubble to first order in fc, 
so it is not completely exact. It can be proven however that, first, the full dependence of this term with k is on the 
form f{k)e'^'^''ki,, being f{k) a function only on the modulus of k, so the structure of the Chern Simons is maintained 
but with a more subtle dependence with the momentum, and second, it was proven by Coleman and Hill that the 
Chern Simons part of the polarization operator given by the one loop contribution (j52p does not get renormalized by 
higher loop contributions, so the coefficient depending on m in (j55p is exact in this model. 

Let us set aside this result for a while and go back to the original problem of graphene in the single cone and low 
energy approximation under the effect of an external magnetic field, but this time adding to the original hamiltonian 
(PU]) a mass term. We can play the same game as we did before, that is, calculate the new Landau level spectrum, 
project the effective action onto this Landau level basis for a finite chemical potential if necessary, calculate the 
average particle number as the derivative of the effective lagrangean with respect the chemical potential, and directly 
read the expression for the Hall conductivity. 

In this case, the Landau level energies are given by the following expression: 



= ± — V2|n|+m2, (56) 

IB 

but this time this relation holds for n ^ only. From the pair of solutions that naively would lead the value n = 
only one of them is actually normalizable (the sign of B determines which solution is normalizable). This means that, 
as in the case of to = 0, there are pairs of Landau levels related by — —En except for the case n = 0, which gives 
an unpaired state. In this case, the average number of particles n is given by the rather cumbersome expression 

" ^ "Sr^^I"'' ~ ^) + ^ (^(M,rn,i3) + ei^, - |m|), (57) 
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where 7 is a integer valued function of the absolute value of the mass, the field, and the chemical potential. What 
is important here is the limit /i -> for the contribution coming from the unpaired Landau level for which the former 
expression reduces to 



n = 

which of course, can be written in a gauge invariant 
ductivity uh'- 



eB m , , 

47r \m\ 

phasion through a Chern Simons lagrangean with Hall con- 



OH = -7--, — r- (59) 
47r |m| 

At this time it can be shown that the gauge invariant current is parity odd (recall the difference with the contribution 
of the zero mode in the massless case: there the magnetic field entered through its absolute value. This is another way 
to see that there is no parity anomaly in t he massle ss limit, although the time reversal symmetry is already broken 
and a genuine Chern Simons term appears) ^ 146l IT50| . 

Comparing the coefficient in front of ([55]) and the value for the Hall conductivity for the lowest Landau level, one 
cannot avoid to wonder if there is any connection between the QHE for massive Dirac fermions and the QED2+1 
theory. 

The connection, if exists at all, must be quite deep, because despite of the superficial resemblance between the 
two Chern Simons theories, some few differences are apparent. In the case of the QHE, both in graphene and in the 
standard two dimensional electron gas, the coefficient in front of the Chern Simons action changes its value when the 
chemical potential eventually crosses a new Landau level, while in QED2+1 the Hall conductivity is fixed. Also, in the 
case of the QHE, the expression for the Chern Simons lagrangean is the ex act on e, while in QED2+1, as we said, the 
momentum dependence in the expression (j55p is valid up to first order in A: [l5ll |. The third and definite difference is 
that in QED2+1 the Chern Simons term is absent if we have pairs of Dirac species related by time reversal symmetry 
(or parity). We can see this in the case of graphene with a mass, where we do have two species of Dirac fermions 
related by time reversal symmetry, so in absence of any time reversal breaking perturbation, the two copies would 
lead to corresponding Chern Simons copies, but with opposite sign, letting a vanishing Hall conductivity. However, in 
the presence of the external magnetic field, both species contribute with the same coefficient to the Hall conductivity, 
so we have a non-zero an- 

We will see that the common root between the systems with a Chern Simons term in their effective ac tion is a 
nontrivial topology coming from the Berry phase. We will also see that in fact playing with the Berry phase (l52j |. we 
can find time reversal symmetric systems with a topological term in their effective actions similar to some extent to 
the Chern Simons term. 



D. The Berry phase and the Integer QHE. 

We have seen in the previous section how we can describe the low energy effective action governing the physics of a 
system (both graphene and the standard two dimensional electron gas) under the effect of an external homogeneous 
magnetic field by means of a Chern Simons term in the lagrangean. We obtained this term under the assumption that 
the system was in the quantum Hall regime, that is, when the magnetic field was so strong to induce the presence of 
Landau levels. The quantization of the Hall conductivity seemed to follow from the existence of such quantized levels. 
Although this picture is true, the presence of Landau levels is not the ultimate reason for the quantized nature of the 
Hall conductivity (although extremely hard to find in condensed matter physics, the example of a single species of 
planar Dirac fermions is the proof of that). Also, in absence of interactions, it is hard to find crystalline solids where 
the time inversion symmetry is absent so we cannot expect to find any term similar to (j55p in any way by looking for 
situations where the one particle spectrum is described by Landau levels, that is, having flat bands in the spectrum 
does not necessarily mean that those bands will display a Chern Simons-like term in the effective lagrangean. 

In order to understand the ultimate reason of the quantization of the Hall conductivity, and envisage any possibility 
of finding similar physics in condensed matter systems, we have to go to the lattice and analyze the same problem, 
but this time not relying on a continuum approximation. We will see that the concept of Berry phase is behind of 
this quantization, that it has a topological nature, and it is the element that eventually will allow us to extend these 
concepts to time reversal invariant systems. 
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1. Berry phase and the quantization of the monopole charge. 

Let us begin this section with a brief review of the concept of Berry phase. We wih fohow a rather pedagogica l way 
and introduce the concept of geometrical phase following the original line of thinking from Berry's original work |152j . 
Berry originally considered the effect on the phase factor of the wavefunction when an slow change was performed in 
the parameters of the problem. 

Consider a quantum one particle system described by the hamiltonian H{s) depending on the set of parameters 
s = (si, .., s„). We will assume that the hamiltonian is time dependent through the parameters s(i). If the parameters 
s change slowly with time, the Schroedinger equation can be solved in the adiabatic limit. If the system of eigenstates 
in principle is not degenerate, at any given t we have 



\n{t)) =e"''^o''-s„(r)+»7.W|„(0)). (60) 

The first term in the exponential is the standard dynamic phase while the second term is a yet undertermined phase 
allowed by the adiabatic theorem. If (l60l) is a solution of the Schoedinger equation H {s{t))\n{t)) = En{t)\n{t)) , then 
the phase factor 7„ adquires the form 



ln{t)^i f^dT{n{T)\^\n{T)), (61) 

or in terms of the parameters s, 

rs{t) rs{t) 

j„{t)^i ds(n(s)|Vs|n(s)) = i / dsAn, (62) 

Js(0) Js{0) 

with the obvious definition An = ("■|Vs|n). The vector An is usually termed the Berry connection. 

The interesting situation comes when we consider the particular case of a system with a multidimensional parameter 
set s and a time evolution where at time t — T, some loop T in the parameter space is performed and therefore 
s(r) = s(0). In this situation, the Berry phase 7„ reads (after using the Stokes theorem) 



i (f> 



dsAn = i / dSVs xAn^i / dS- Bn- (63) 
Js Js 

In the absence of any singularity in the parameter space or in the vector field ^„(s), any loop T can be continuously 
deformed to zero and there is no Berry phase effect in the wavefunction. The funny situation comes when the 
parameter space is itself nontrivial or the field An shows some singularities. Let us recall that although we have made 
use of the adiabatic hypothesis, it is actually not needed since the time t (or r) disappears from the expressions (|6ip . 
(|62p. and (|63p and all the expressions are written in terms of geometrical quantities defined in the parameter space. 

We can follow the original paper of Berry and transform (1631) a little bit more: 



ln = i j rfSVs X (n|Vsn) = i /dS(Vsn| x [Vsn) (64) 

= i/dS^„_^„(Vsn|TO) X (m|Vsn). 

If we consider the situation where the eigenvalues En are not degenerate, for n ^ m, we have (£'„ — E„-i){m\\I s\n) = 
{m\V sE[\n) we find the important final relation; 

7„ = -^ / dS 2^ — —2 . (65) 



s 



, {En — EmY 

n^ra 



Although in the expression (1651) we have written the Berry phase 7„ making use of the non degeneracy of the eigen- 
states of the system, the previous expressions can be used to calculate this phase in the situations when degeneracies 
in the parameter space exist. We will use the low energy hamiltonian in graphene to illustrate this situation and as a 
toy model to show that the Dirac quantization of the monopole charge is actually a topological quantization due to 
the existence of a nontrivial Berry phase in that system. 
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This time we will take into account both species of massless Dirac fermions in the low energy sector of the spectrum 
in graphene (our parameter space here will be the momentum space). As usual, the effective hamiltonian around the 
Fermi points if, K' can be written as 

Hr^k' = V [r^axK + IrCyky) . (66) 
The eigenstates around both Fermi points are E{k) = sv\k.\, and the eigenstates are collectively written as 



where s = ±1 and t — ±1 stand for band index and Fermi points, respectively, and 0k is the angle defined by the 
wave vector k with the horizontal axis. The Berry connection ([62]) is easily calculated to be = 1, 2) 



In terms of the Berry curvature, we get 



^lr{^) = V-W- ^^^^ 



Bl^ = ^S(^Hk), (69) 



that is, both species carry a quantized monopole charge of value one half (sitting at the K and K' points), and it has 
opposite value for each specie. In this particular case, the total magnetic charge i s the sum of the two charges in (p^ 
and this zero, being the Nielsen-Ninomiya theorem behind this particular result |l53j |. We have to make a comment 
here concerning the particular value the monopole charge takes. It is a general statement that around degenerate 
points the monopole charge takes half integer values. We will see in the next sections that in general this topological 
charge takes integer values instead. There is no contradiction since, as we have said, the Nielsen Ninomiya theorem 
ensures to have an even number of degenerate points, and the total charge will be zero, if the degeneracy points are 
related by parity or time reversal invariance, or some integer, because the sum of even number of half integer charges 
is an integer, if the points are not related by theses symmetries. It is common to find these degeneracies lifted by 
some allowed perturbation. When the total charge is zero or an integer depends on the nature of this perturbation. 

This case exemplifies how a nontrivial structure in the wavefunctions leads to a vector field (Berry curvature) which 
is singular at some points of the parameter space, and the total integral of the curvature associated to this singular 
connection over a closed surface (in this case, any two disjoint spheres containing separately the points K and K') 
leads to quantized value. 

2. QUE on the lattice and the magnetic Bnllomn zone. 

Let us consider, for this time, fermions defined on the lattice whose dynamics under the effect of a magnetic field 
are described by the following tight binding hamiltonian with coupling to nearest neighbours. In what concerns the 
basics, ther e is no much more complication if we study the prob lem in the square lattice or in the honeycomb lattice 
instead jl54l |. We will follow the original work of Hofstadter [l55j |: 

H^Y. t{^^.^■J)C+{^^)C{Ii■J), (70) 

<i,j> 

where Ri labels the sites of the square lattice, and the four nearest neighbours are Rj = {R^ ± aea;,Ri ± aey\. The 
hopping term t(Ri, Rj) is defined by the standard Peierls substitution: 

t(R„R,) =te*^^>^'^^'^^''''. (71) 

In the Landau gauge A — i?(0, a;, 0), the phase in (|7ip is just eBax, where x = am labels the a;— component of the 
position Rm and a is the lattice spacing. Due to translational invariance along the y direction, we can rewrite the 
hamiltonian H now reads 
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^ = i Ek, C+{ky)Cra+l{ky) + C+{ky)C„,^l{ky) + (72) 

2 COS {2iT4>m + kyo) C^{ky)Cm{ky), 

where (p is the magnetic flux within each plaquette. When this magnetic flux per plaquette is a rational number, 
= ^, the cosine term in (1731) is invariant under the change m — > m + (7 and so does the hamiltonian. The Bloch 

theorem thus tell us that Cm+q = 6*'^'="'?'' Cm for q fixed by the magnetic flux. Our Schroedinger equation is a q x q 
matrix problem, with q eigenvalues E = -E'j(k) with j = 1, ...,5. 

We can fix the periodicity of Ej (k) by noting that the points k^ and + ^ are equivalent in ([75)) so Q < k^ < 

< fcy < so we can define a Brillouin zone, which has the form of a torus (that the Brillouin zone has is a 
torus and thus a closed surface is one of the key ingredients for the quantization of the Hall conductivity, as we will 
see) . The important message here is that although the presence of an external magnetic field prevents the system to 
enjoy the translational symmetry of the original system at i? = 0, we still have an enlarged translational symmetry, 
which is the same as to say that the original Brillouin zone is split up in q copies. Also we have enlarged the number 
of bands and the new eigenstates are eigenvectors of q components. 

Instead of studying the consequences of this effective folding by the general properties of the solutions of ([751) , for 
our purposes it is enough to analyse two specific simple examples for two values of the pair (p, q) . 

Let us start for the pair {p,q) — (1,2). In this case, the Schoredinger equation is a 2 x 2 matrix equation, which 
takes the particularly simple form of 



H = 2t(-^''<fy} cos(aM\ 

\ COS(aK:r) COS(aKj^) J ^ ' 

It is worth to mention that the situation with half of the flux per plaquette is known as the " pi-flux phase" . It is also 
worth to mention that for this particular choice of the flux the system is actually time reversal invariant, because as we 
said before, there is a periodicity in (f) and the system with = i is equivalent to the system with </> = — i. Nevertheless, 

the example g = 2 is interesting because the spectrum takes the form i5±(k) = ±21 ^ cos{akxY + cos{akyY. 

The spectrum consists in two bands that are degenerate at the points (7r/2, ±7r/2). Around these two points, the 
hamiltonian (|73p takes the familiar form: 



H = vTz(Tx6kj; + vazSky, (74) 

that is, the low energy hamiltonian of the system consists in two species of massless Dirac fermions, which are, 
interestingly, related by parity. This example illustrates a general property of ([75)) for even values of p: the band 
structure consists in gapped bands except for the two central bands that are degenerate at q isolate points in the 
Brillouin zone (as in graphene, labels the t wo D irac species). Around those points, the low energy system can be 
described by q species of massless Dirac Doints (l56j |. For odd values of q the bands are globally separated by gaps in 
the entire Brillouin zone. In the case for {p,q) = (1,3) the hamiltonian is the following 3x3 matrix 

/ cos{aky + 27r/3) 1 e^**^-" \ 

H = t\ 1 cos{aky + 4:n/3) 1 (75) 

\ e-^'''-" 1 cos{aky + 67r/3) / 

The eigenvalues i?(k) are the solutions of the secular equation 

E'^ - &t^E = 2t^{cos{3k^) + cos(3fcy)) 

and they are depicted in fig. |6] (right) |l56l Il57l |. 

These two examples illustrate to what happens when electrons living on a lattice experience the effects of an external 
homogeneous magnetic field. The original Brillouin zone is split up in q copies of a new magnetic Brillouin zone, 
which is still a torus. Also, the new wavefunctions become multicomponent Bloch functions juk)- This fact seems to 
be rather paradoxical if we again compare with the problem of Landau levels. In the previous sections we found that 
the spectrum of electrons in two dimensions under the effect of a external magnetic field consisted in discrete Landau 
levels, and, due to the explicit dependence of the vector potential A with the spatial coordinates, the invariance 
under translations was lost. A deeper look on the problem will tell us that this is not the case. Effectively, the 
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FIG. 6: Bandstructure of the pi-flux phase corresponding to the value q = 2 (left) and q — S (right). 

physical momenta 11 satis fy [na;,!!^] = —ieB, so two elements of the translation group do not commute in general. 
If T(R) = e™, then[l5| 

r(Ri)T(R2) = T(R2)T(Ri)e*'=^-(^i^^^). (76) 

They do commute, however, when the phase $ = eB • (Ri x R2) is an integer number, that is, when the flux passing 
through the area defined by Ri x R2 is an integer. In terms of the original primitive lattice vectors, this condition 
means that if the vectors are written in the original lattice basis, Ri = pa and R2 = gb, we have the condition for 
flux per unit cell to be a rational number, pq(j)piaquette = n, being n some integer. This condition is precisely the one 
that we found at the beginning of this section, when we considered the eigenstates of (|73l) . We then find that even 
when the original translational invariance is lost, we can still label the eigenstates within each Landau level in terms 
of a wave vector k when the magnetic flux per unit cell is a rational number. 
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3. Berry phase, Chern number, and the Hall conductivity. 



In this section we will close the circle and see the explicit connection between the Hall conductivity and the Berry 
phase constructed through the Bloch states corresponding to the new magnetic subbands Ej{k). Our task is to find 
an expression that relates the Hall conductivity and the Berry phase 7„, and second, we need to demonstrate that 
the Berry connection defined by the Bloch eigenstates of (l73l) is a singular vector field at some points of the Brillouin 
zone, and thus the total monopole charge enclosed by this surface is an integer. 

Let us employ linear response theory to find the desired relation between an and t^ IISSI ]. Assume that, together 
with the external homogeneous magnetic field B applied perpendicular to the system, we apply a weak electric field 
E pointing, say, in the x direction. The corresponding vector field thus takes the form = E^x and A ~ (0, Bx, 0). 
We shall assume for simplicity that the magnetic subbands are globally separated by a gap in the entire Brillouin 
zone, so Ei{k.) ^ Ej{k.) for every momentum k and i ^ j- Using the magnetic Bloch functions in the real space 
representation |V'n(r)) = e*'^""|u„(k)) as the basis of unperturbed wavefunctions, simple perturbation theory gives the 
first correction for the perturbed wavefunctions: 



|.„(k))(^) = |.„(k)) + ^eE. ^(,),^(,) • (^^) 

n^ra 

In one particle band hamiltonians the velocity operator is defined as the gradient of the Hamiltonian, v(k) = 
Vki?(k) and the averaged electronic current is just the averaged velocity times the electric charge: (j) = e(v). We 
then can use the expression for v(k) and the perturbed eigenstates ([771) to calculate the first order correction to the 
y component of the electronic current of the state n (as long as (v„) is a real function, we must include also the term 
given by the complex conjugate of ([77])): 



\vy\un) + leE^ j ^ _^ + c.c. (78) 



Again, we can make use here of the identity (En — Ej^^){m\V\i_\n) — {Tn\V\i_H\n} as we did two sections ago and 
plug it in ([75|) giving 



^ {Um\dk^Un){Um\dkyUn) - (Um | ^fc^ Wn) (Urn 1 9fc, U„) . (79) 

We readily read in (j79p a term similar to th e on e in the last part of the expression so the second term of ([7^ 
can be written in terms of a Berry connection (l52j | 



If the Fermi level lies on a gap above the band -E„(k), the first term of the right hand side of (|80l) do not contribute 
to the total current. Because we have the relation (j'*) — e(u*) and (j*) = cTijE^ , we read from the value of the 
Hall conductivity: 



e2 ^ 



an = rr.y - — 2^ 



d^kiVk X Anik) = £3 j( yl„(k). 



where T is the loop defined by the corners of the magnetic Brillouin zone. This is the desired relation to prove: the 
Hall contribution to the conductivity is the total fiux passing through the Brillouin zone of a Berry curvature. When 
we discussed the quantization of the charge of a magnetic monopole in the context of graphene, we found that if the 
Berry connection was singular at some points inside some closed surface, the monopole charge was quantized. Here 
we must show that the Berry connection in (|8T|) defined through the magnetic Bloch eigenstates |V'n(r)) is singular 
at some point within the Brillouin zone. 

The proof that there are points on the magnetic Brillouin zone were th e Be rry connection is singular is rather 
involved and we will only sketch it, referring the original work for details [33l Il59j |. 
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We have stressed in the previous section that the magnetic Bloch wavefunctions that are solutions of ([73)1 are no 
longer eigenstates of the translations of the original lattice, but they are of the new enlarged unit cell, that contains 
an integer value for the magnetic flux passing through the new magnetic plaquette. If we consider a loop consisting 
in such of one of these magnetic plaquettes, the phase of the magnetic Bloch wavefunction |V'k(r)) = pk(r)e'^'''^''-' 
will have changed 27rp: J (ilVi6'k(r) — iirp . This is nothing but the Aharonov-Bohm effect of an state encircling a 
region of the space where there is a nonzero magnetic flux. It means that the wavefunction has the form of a vortex 
at some point of the magnetic unit cell in order to accommodate the nonzero phase in a closed loop. As long as 
the wavefunction is not singular on its arguments, the wavefunction must be zero at the point where the vortex is, 
and for some point in the magnetic Brillouin zone, as long as the phase 0k(r) is a function of k. Because the Bloch 
wavefunction has zeros on the torus defined by the magnetic Brillouin zone, the Berry phase will have points where 
it is ill-defined. 

Let us now assume that we can divide the magnetic Brillouin zone in two regions Rl and R2 separated by any loop 
7, one containing the zero of the wavefunction. In the region that no contains the zero, the phase of the wavefunction 
is well defined globally, but in the region that contains the zero the phase might pick up any arbitrary value. Once 
chosen, we can allways choose the phase continuously in a small neigborhood around the zero with the condition that 
the phase cannot be defined globally. When we compare the wavefunctions away from the zero, the two choices of 
phase must be consistent through a gauge change, that is 4'/ — e^^^^^'i'ii, or, equivalently 

Anj(k) = AnjiCk) + Vk/(k). (82) 

By using the Stokes theorem, we can apply the definition of Hall conductivity (j8ip and write each contribution for 
the two regions: 

2 r 2 



<^H = —^ d^kVk X An.i + / d^kVk x Anji, (83) 
47r Jjii 47r J_fj2 



and then, because Rl and R2 have opposite orientation, 

<^H = ^ J^iAnJ - AnJl)dk ^^J^ Vk/(k)dk. (84) 

The expression (|84p means that the Hall conductivity counts the number of times the function /(k) winds around 
27T as k travels along the loop 7. Since this winding number does not change when the loop 7 is deformed (and also 
it contains the zero of the wavefunction), this winding number is a topological invariant, known as the first Chern 
number, and it takes integer values, which is the result we wanted to prove. We have written the Hall conductivity 
as the surface integral of a Berry curvature, which is nothing but the total Berry phase calculated along a loop 
corresponding to the "boundary" of the Brillouin zone. In order to have a nonzero an, the wavefunction |'0„(k)) 
must pick up a nonzero Berry phase when circulating along the mentioned loop. If the wavefunctions have no Berry 
phase, or equivalently, the Stokes theorem applied to this situation tell us that if An is not singular and we perform 
the integration over a surface without boundaries, the integral is zero. This is what happens with most of the band 
insulators. However, there is an intermediate case, which is when the system has an even number of singular points 
in the Berry connection, by requirements of time reversal symmetry, and also the winding number of those points 
have opposite value (we saw this in the case of the low energy hamiltonian of graphene). The total Berry phase is 
zero (and so does an) but we are still able to find a modified topological quantity which takes into account the value 
of each singular point separately. We will analyse this situation in the next section. 

Another interesting observation concerning (j81|) is that when in normal band insulators the Fermi level lies on 
the gap, the Bloch wavefunction of the groundstate is made of linear combination of (maximall y) loc alized Wannier 
functions, that is, the wavefunction of a state at the point decays as (or faster than) |r — r, |~^ [l6C| . It means that 
the Bloch wavefunction |'0(k)) = e''"'|u(k)) in the momentum space is an analytic function of k. A Hall nonvanishing 
Hall conductivity means that we cannot define globally the phase of the Bloch wavefunction, so it will not be analytic 
[l6ll |. It thus implies that the ground state of a Hall system cannot be described in terms of maximally localized 
Wannier functions. This observation has interesting consequences in the mode rn ge ometrical theory of the electric 
polarization based on the Berry phase of maximally localized Wannier functions |162|. 



E. Topological insulating phases in graphene. 



In the previous sections we have found a relation between the Hall conductivity an and the Berry phase for the 
Quantum Hall effect. We have also seen that despite of the presence of an homogeneous magnetic field, we can still 
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define a Hamiltonian which enjoys translational invariance (this translational invariance is not the same as the original 
system) and with the magnetic Bloch wavefunctions we can define a Berry connection, whose circulation around the 
perimeter of the Brillouin zone gives the quantized value of <jh- However, at the beginning of the previous section, 
we compared the QHE with a single specie of massive Dirac fermion in two spatial dimensions finding in the two 
cases a Chern Simons term in the effective electromagnetic field theory, being the coefficient of the CS term the Hall 
conductivity. Let us see how we get the same conclusion for QED {^2+1) by using the Berry phase ((65|). We will make 
extensive use of this calculation in the subsequent sections. 

The hamiltonian corresponding to the action (j5ip takes the following form, in the momentum space {v is the Fermi 
velocity and cr = [dx, o'y))'- 

H ^ verk + maz. (85) 



The eigenvalues of this hamiltonian are, of course E+ = —E- = yV|kp+m^ = _Ek, with the corresponding 

eigenvectors |+) \[^^^ (l, iK^Ezf ^^'^ = ^^3^ {-1^^^ l)^- We have used the notation k = + it 
and k* = kx — iky. Because we have two bands, the Berry curvature takes the simple form, for the positive band: 



y 



= * -^2 ■ (86) 



After pretty simple algebra, we get 



Bz,^^-"- ^ (87) 

2 (w2|k|2 + m2)2 

From (|86|) it is easy to find Bz,+ — —Bz-- We then end up with the expected result ([55]) (as usual fi — 1): 

2 J 47r^ (w2|k|2 + ^2)2 47r |to| 

We have closed the circle and demonstrated that the common origin of the Hall conductance in the case of the 
QHE and QED2+1 is the nontrivial Berry phase in the momentum space picked by the eigenstates of the system. 
To end this introduction, we can comment on the properties of y8z(k) under the effect of discrete symmetries, like 
time reversal symmetry and parity. If the system is time reversal invariant, one can show that Bz(—i^) = —BzCk) 
while if the system is parity invariant the Berry curvature satisfies (—k) = Bz(k-)- These two relations imply 
that if the system is time reversal invariant, the Berry connection does not need to be zero, but the integral of the 
Berry curvature, i. c., the Hall conductivity, vanishes. Also, if the system is both time reversal invariant and parity 
invariant, — ^^(k) = Sz(— k) = S^(k), so the Berry curvature is zero. 



1. Haldane model and the valley Hall effect. 

Now it is turn to look for situations different to the already analized QHE where the system displays similar physics. 
The condition of vanishing Hall conductivity for time reversal invariant systems might seem too strong and one might 
not expect the physics described by the Chern Simons term in any time invariant system. However, with the aid of 
graphene as an ideal laboratory to play with, we will see that this is actually not the case and we can still find time 
reversal invariant system displaying nontrivial topological f eatures. F or the details and a for an exhaustive list of 
references we refer to the excelent reviews on the sub i ect [lOl . 1 1 63l . . 

Let us consider how we can modify the lattice hamiltonian describing the low energy spectrum in graphene. We 
know that, if we do not consider the real spin, the low energy hamiltonian of graphene is described by the hamiltonian 
(l66l) . If we insist in opening gaps for the whole set of bands, we can add to (|66|) the following matrices. Mi = l-rCz, 
M2 = Tzdz, and the two matrices M^^^ = Tx^yCx- Physically, Mi represents an atomic energy invalance between the 
atoms belonging to the two sublattices in the honeycomb lattice. In graphene it might appear due to a mismatch 
between the graphene sheet and the substrate lattice or in the boron nitride. M^^^ physicaly correspond to two 
degenerate realizations of the coupling between the electrons and a Kekule distortion (a kind of optical phonon in 
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the honeycomb lattice). The Kekule distortion is a quite interesting coupling by its own right. In its vector form 
Hjiekuie = •^(r)TCTa; it constitutcs a vector coupling with the electronic structure at low energies. This coupling iis 
symmetric under chiral rotations of the form U{6) = e"'^^"^ if the A(r) rotates as a vector in two dimensions. If 
the coupling Hxekuie is generated by any form of interaction, the particular choice of A(r) implies an spontaneous 
chiral symmetry breaking in graphene. As long as we are at finite temperature, increasing T would eventually lead 
to a chiral symmetry restoring by the appearance of defects in the vector field A(r), through a Kosterlitz-Thouless 
type of transition, which has some topological concepts behind, as we commented in the introduction. In terms of 
the Berry phase we are discussing here, this coupling is trivial. 

The first, third and fourth matrices are time reversal invariant, as long as the time reversal operation acts in the 
following way on the tight binding hamiltonian: &H{'K)Q~^ = K), that is, it interchanges the two Dirac species 

and takes the complex conjugate. The only matrix that breaks time reversal symmetry is Haldane devised a tight 
binding hamiltonian which in the low energy sector a mass like M2 appears, by adding complex second neighbour 
hoppings which opposite value in the two sublattices and an unbalance site energy between the two sublattices of the 
honeycomb lattice. 

Let us make extensive use of the calculation performed in the previous section. The haldane model in the low 
energy sector reads: 



H = VTz<7xkx + vlrCTyky + TOT^CT^. (89) 

Because there is no mixing between the two Dirac species, we can apply the formula (j86p to each specie separately. 
For the first specie we find (|87p . The important part comes now. The eigenstates for the two species are related 
by |2,+) = iay\l,+) and |2, — ) = ~iay\l,~) while the velocity operators satisfy wj'^ = ayV k.H'^'^(}i)ay, where 
the superscript labels the two species. Taking into account the relations between the eigenstates and the velocity 
operators, one finds S^^^(k) ~ S^^|^(k). It means that after integrating over the momentum domain, both species 

equally contribute to the Hall conductivity, and we get uh — ^sign{m). This is consequent with the fact that the 
system is not time reversal invariant. 

What if we consider now the same model, but with Mi instead of A/2? This time the eigenstates for the two species 
satisfy |2,+) = (Ta;|l,— ) and |2,— ) — cr2;|l,+), that is, this time the positive and negative solutions for one specie 
are related to the negative and positive solutions for the other specie, respectively. Now, in terms of CTa,, the velocity 
operators for the two species now satisfy a^Vjll^'^'^ax = —V jH'^^\ so, putting all together in the expression for the 

Berry curvature, we find B''^\_ = — that is, now both species have opposite Bz,+, so the total Berry curvature is 
zero. 

It is interesting to note that, despite of the total Berry curvature and uh are zero, the contribution of each specie to 
Bz.+ does not. So, the two species have a nontrivial topological structure separately represented by a nonzero Berry 
curvature, as long as there is no coupling between them. Then the question is if there is any physical observable 
depending on the separate Berry phase for each specie. In order to ask this question, we will make use of the 
Chern-Simons effective action formalism. 

Let us denote the Hall conductivity of each specie by gk and uk' ■ We know that uk+'^K' = so uk' = — cta' = ^o". 
Let us postulate the existence of two independent U{1) gauge fields, one being the usual electromagnetic field Ai and 
a second one Vi. The gauge field Ai will be coupled to the standard fermionic current J'^, while the vector Vi will 
couple to a quiral current J|', that couples with opposite sign to each Dirac specie. The gauged Dirac lagrangean now 
reads: 



£ = (7°^ - «7k - ml) Vk - e J^A^ - eJt^V^,, (90) 

with the obvious definition of the 7-matrices and ij^ = ip'^j^ = ip'^ Mi. We can play the game done previously in the 
context of the QED(2+i) and calculate the effective one-loop action for the two U{\) gauge fields integrating out the 
fermion fields.. Among all the terms in the effective action, there are two nonvanishing terms that have the structure 
of (pseudo) Chern-Simons terms: 



Spcs - ^{(TK -<JK')J rf^xe^'^" {V^d.Ap + A^d^Vp) + eJ^'Ap + eJ^^^, (91) 

where the couplings with the fermion currents are added by hand. The average of the quiral current is the found by 
taking the functional derivative of (|9T|) with respect the chiral field V^: 
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W) ^ {J'k - J'k') = ^ = ^oe^-^'d^A,. (92) 

This expression means that there is a nonzero transverse current reacting to an external electromagnetic field, which 
is actually an inbalance between currents belonging to different Dirac species. In order to make this statement more 
precise, let us employ the Callan-Harvey mechanism discussed above. Let us take linear combinations of the fields 
Af^ and V^: = A^ + and = — V^, so = 5(0^ + and V^, = \{a^ — &^), and we do the same for the 
currents: = + J^, and J!f = J'^ — J^,. Writing (j9T|) in terms of and we have: 



Spcs = j d^xe'"''' {af^d„ap ~ b^di,bp) + eJ^a^ + eJ'^,bf^, (93) 

That is, two independent Chern-Simons terms with opposite sign. We note that the corresponding Hall conductivity 
is the Berry phase for each single Dirac specie. Now put the system described by (|93p in a finite region. Applying 
the Callan-Harvey mechanism, we find that the system develops two counterpropagating chiral currents belonging to 
different Dirac species. This is known in the context of graphene as the valley Hall effect. 

In the case of the QHE, we found that the system developed chiral one dimensional propagating states at the 
boundaries in a number given by the value of ajj. Because the Hall conductivity can have any integer value, we can 
have the situation of getting any integer value of chiral states. In the case of QED^2+i) and the valley Hall effect, 
this is not possible, because ajj can only take two values <th = ±1 in appropriate units. This is another important 
difference between the QHE and QED(^2+i) and time reversal invariant topological insulators: the topological invariant 
in the two later cases is a Z2 invariant. We note by passing that in the Haldane model, the topological invariant is 
also a Z2 invariant, and we have one chiral state propagating along the same direction at the boundaries. 

We can make a comment concerning the robustness of the edge states in these systems. In the case of the QHE, all 
the edge states present in the system are propagating in the same direction, independently of their number. It means 
that those states cannot scatter backwards, so they are dissipationless and the Hall conductivity remains quantized 
even in disordered realistic samples. The same happens in the Haldane model. In the case of the valley Hall effect, 
however, we have pairs of counterpropagating edge states, but belonging to different Dirac species, so in principle 
the edge states can backscatter if there is any perturbation or interaction coupling the two Dirac species, so the 
system is not protected against all kinds of perturbations, as it happens in the QHE. However, it turns out to be that 
the backscattering rate between species is extremely small as it can be seen in the case of impurities (cite paper of 
paco/gladys) and interactions (cite me/Henning). We then conclude that we can find out the way of unraveling the 
nontrivial topological structure of some time reversal invariants by using the Berry phase. Here we took the route of 
analysing the system in terms of an effective pseudo Chern Simons theory, but, as it happened with the QHE, we can 
study the topological structure of the system by means of the behavior of the Bloch wave functions in the Brillouin 
zone. This has been done in (cite all). We will not pursue this line here. 



2. Quantum spin Hall effect. 

In graphene, when spin is taken into account, the possibilities of mass-like couplings notably increase. Among 
all of these possible couplings, there is one that is of paramount importance, since is the on that will be generally 
responsible of the Quantum Hall effect in general time reversal invariant materials, different from graphene. The 
physical effect involving spin that will lead to the coupling we are talking about is the spin-orbit coupling. The 
spin-orbit coupling, as its name points out, couples the spin momentum with the orbital motion of the electrons in 
the lattice. It takes the general form of Hgo = AS • L, with A some constant depending on the microscopic details of 
the atomic potential. It is easy to see that the term Hso respects time reversal symmetry because both operators S 
and L flip their directions under reversing the time. The low energy effective spin-orbit coupling in graphene takes 
the form (we refer the literature for explicit lattice derivations: cite dresselhaus(1965) and others): 



Hso = AsoCTzTzSz- (94) 

It is interesting to note that the matrix structure of (|94)) is similar to M2 in the Haldane model. In the Haldane 
model, the mass does not respect time reversal symmetry, but, because the spin changes its sign under time reversal 
invariant, the product M2SZ = cr^TzSz as a whole is time reversal invariant. It is then clear that the low energy 
hamiltonian for graphene in addition to the term (j94[) consist in two time reversal related copies of the Haldane model. 
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From the discussion of the Haldane model, we infer that the term ([94]) opens a gap for the two spin components, and, 
because the effective field theory in the case of the Haldane model possesses a genuine Z2 Chern Simons term, the 
low energy effective field theory in the case of graphene with an spin-orbit term will have a term similar to (|93|) , but 
this time each current does not correspond to different Dirac species, but they belong to different spin projection. It 
means that in a confined geometry, the spin-orbit term will induce at each boundary, two counterpropagating edge 
states with opposite spin projection, so if one has a wide sample and is capable to manipulate the system at the 
scale of the edges, under the effect of an appropriate bias voltage, he or she will lead an spin filtered current, which 
constitutes the holy grial in the field of spintronics. 

Despite of this so much appealing property, the spin-orbit constant Ago in graphene turns out to be extremely 
small. In units of temperature, A^o is estimated to be of the order of 0.01 K, so any other effect in graphene is likely 
to have a bigger characteristic scale and be dominant compared with -ffs-o- Also, in graphene there is other coupling 
that mixes spin and orbital movement. It is known as the extrinsic or Rashba spin-orbit coupling, which has the form 
Hr = Xr{s X p) ■ ez- The Rashba coupling takes the following form in the low energy Hamiltonian of graphene: 



Hr = Xr {r^axSy - IrO-ySx) ■ (95) 

The Rashba term, depending on and Sy competes with the spin-orbit term (j94p . and it can be found that the 
system remains to be an insulator ii Xr, < Ago, so while this condition is fulfilled, the spin Hall conductivity remains 
nonzero and the system is a topological insulator. We have to m entio n that in this case, although the spin Hall 
conductivity is nonzero, it does not generally take quantized values jl65j |: 



An Xr 



Xr + As 



Xr 



(96) 



It is interesting to note that in order to induce backscattering between the two spin filtered edge states at the 
boundary, we must have a perturbation that couples the two spin polarizations, that is, a time reversal breaking 
perturbation. This is a manifestation of the fact that the edge states are protected against backscattering by time 
reversal invariance. 

The hamiltonian (I94p together with the continuum effective Dirac hamiltonian for graphene are useful for describing 
the physics we are interested to in a simple way, but of course not all the two dimensional topological insulators must 
be described at the lattice level by this hamiltonian, although the effective field theory contains a Chern Simons 
term. We saw an example of this in the lattice description of the QHE. Indeed in the case of graphene, the discussion 
made so far remains valid as long as the value of the gap is not too big compared with the bandwidth, and the 
Berry curvature concentrates around the Dirac points. For these situations where the Dirac approximation is not the 
correct one, one may take to routes: to construct a fairly good effective hamiltonian for the electrons that captures 
the topological features of the system, or to develop a general theory for the Z2 topological invariant in time reversal 
invariant systems taking into account all the information contained in the Brilloin zone. 

The second route is being subject of current research and several methods have been prop osed to construct such 
topological invariant. We will not review them here and we refer to the Iiterature |l66l - [l70 |. Concerning the first 
route, we will use bilayer graphene as an example. 



3. Topological insulating phases in bilayer graphene. 

For bilayer graphene in the Bernal stacking the unit cell is made of two pairs of carbon atoms, each pair belonging 
to a different graphene layer. The simplest lattice hamiltonian that shows nontrivial topological features is the 
one where only hopping amplitudes are allowed between the most adjacent atoms in the two layers. If we define the 
wavefunction of the system for the low energy effective model around the K point by — (ai.k, &i,k, 0'2,k, ^2,k)"^ being 
the components the amplitudes in the sublattices a and b of the two layers, this simplest model reads H = ^k'^k^k 
with 

where we have used the standard definition, a± ~ {ax Hcy) /2, and t± is the hopping amplitude between the 
sublattices 02 and bi. The effective model for the other K point is found from (|97p by substituting by —k^- In 
what follows we will use the following convention for the different Pauli matrices appearing in the model. The cr 
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matrices stand for sublattice, r stand for K' points, /x stand for layer index, and s refers to the spin quantum 
number. The system consists then in four bands around each K point and are degenerate in spin, i.e., sixteen bands 
for the whole effective low energy model. 

In absence of any other coupling or perturbation and for energies much smaller than t±^ the model can be simplified 
further to an effective two band model for each K point. As described in Sect. Ill Bl this effective model for the 
amplitudes (02, k, ^i,k)"^ is 

n^! = j- (98) 

In some cases the model (|98p is enough to analyze the topological structure of the system. This is the case, for 
instance, when an external homogeneous magnetic field is applied. The low energy part of the Landau spectrum 
of the system can be derived directly from ([M)) to get En = ±0Ji\/n{n — I) [53], that is, each valley and spin 
projection contributes with two zero modes instead of one, as it happened for monolayer graphene. This change can 
be experimentally seen in the value of the Hall conductivity an for bilayer graphene at zero carrier concentration, 
being the value of (Jh twice the value in monolayer graphene |171j . 

Another interesting effect is the application of electrostatic potential difference between t he tw o layers. This effect 
is modeled by the term V = VoIsMzItIo- in ([97]) ot: V = Voaz in the effective model ([98|1 ^172|. By making use of 
([62|l and the eigenstates of T^k + it can be shown that the Hall conductivity of bilayer graphene in presence of a 
voltage difference is zero as long as this voltage difference respects the time reversal inversion symmetry, but as it 
happened with monolayer graphene, we can still have a valley Hall effect with value of the valley Hall conductivity 
to be aK = -^sign{Vo). In this case, it means that after applying a voltage difference between the layers a pair of 
counterpropagating edge states polarized in valley appear at the edge of the sample. This is quite appealing because 
is by far easier for experimentalists to apply an external voltage difference between the layers in bilayer graphene 
than opening a gap of the type Mi in monolayer graphene. Concerning spin, we expect to find similar topological 
structure than in monolayer graphene when the intrinsic and Rashba spin orbi t term s are considered now. The most 
general term involving both effects in bilayer graphene can be found to be [l05l . flTSt : 



Hso Xllf^TzS^a^ + \2iJ-zTzSz^cr + A3 {pizlrSxiyy - fJ-zTzSyCTx) + A4 {iMylTSxCTz ~ fJ-xTzSy(Tx) ■ (99) 

In the hamiltonian (|99p the term with Ai is the equivalent version of the spin-orbit term in monolayer graphene, 
while the term with A2 is proper from bilayer graphene. Both terms open nontrivial gaps in the bilayer spectrum 
and lead to a quantized value for the spin Hall coefficient, as — ^sign{Xi + A2). The terms associated to A3 and 
A4 are Rashba terms which destroy the quantization of the spin Hall response in a similar manner that happens in 
monolayer graphene. 

An interesting situation appears when the intrinsic spin-orbit terms and the V term are considered simultaneously. 
The system displays both Spin and Valley Hall effects, being both Hall conductances quantized, 

as ^ — {sign{Xi + A2 + Vq) + sign{Xi + A2 - Vq)) , 
Ztt 

and 

aK = TT (si5"(Ai + A2 + Vo) + sign{Vo - Ai - A2)) . 
Zn 

As usually happens, the quantized values for the distinct Hall responses are associated to the fact that a gap is opened 
in the system. Of course bilayer graphene is not an exception and the presence of a finite Fermi surface makes any 
Hall response to be nonquantized. 

To finish this section, we can comment on the fate of the edge states when gated graphene bilayer is considered. As 
we have said the system shows a nonzero valley Hall conductivity ax so if we apply the Callan-Harvey mechanism 
to this situation we will fi nd two conuterpropagating edge states (without taking into account the real spin) filtered 
in the valley number |l74| . As it happened in monolayer graphene or in the QHE, the topological number gives us 
only the number of propagating quiral states, but not their velocities, being them dependent on the microscopic 
details. It might well happen that both counterpropagating chiral states posses different velocities. This fact have 
important consequences for the transport properties, because the backscattering rates when disorder or interactions 
are considered strongly depends on these values |l75| |. 
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VI. SUMMARY AND FUTURE DIRECTIONS. 



One of the purposes of this review, besides presenting some QFT aspects of condensed matter in a language available 
to high energy colleagues, was to explore the theoretical aspects of the problem that can still be developed in the near 
future. The giant explosion of graphene publications after its synthesis has explored all possible aspects of it and the 
latest activity is mostly centered on the applications, devices, and exploring standard condensed matter phenomena 
that are no very different in graphene than in a usual electron gas. 

It is fair to say that most of the special properties of the material originate at the Fermi points and as such are 
difficult to access. New suspended samples allow for the measurement of the electronic properties within a range of 
p « 10®cm~^ 01 ep ~ ImeV in single layer graphene. The density dependence of the effective mass and Fermi velocity 
at low carrier concentrations showed an enhancement of the Fermi velocity of almost a factor of three. This effect 
is consistent with the expected renormalization of the Fermi velocity in the absence of screening, given a graphene's 
"fine structure constant" a = /vp ~ 2.3 — 2.5. Similar estimates for a have been used to explain the observation of 
plasmaron satellites in samples with high carrier concentrations. This value of a puts single layer graphene into the 
intermediate to strong coupling regime, although its flow towards zero takes the system to the weak coupling regime 
at low energies. 

The issue of opening a gap in monolayer graphene is one of the few fundamental questions that still remain 
controversial in the topic. Some recent lattice gauge theory calculations indicate that the physical parameters of 
graphene can be at the edge of the symmetry broken region. Given the interest of the subject on the theoretical 
side and even more for the practical applications, a special effort has to be devoted to develop more sophisticated 
calculations in graphene. We must although warn that there are at the present no experimental evidences favoring a 
gap opening from interactions in the intrinsic monolayer system that has been tested at incredibly low energies 85] . 

There are some QFT aspects of graphene of recent interest that have been left aside of this review not because they 
are not interesting and worth to pursue but for lack of spacetime, or of expertise, or because we consider that the field 
is not yet mature. A very interesting aspect that is being obj ect o f great attention and will probably be developed in 
th e nea r future is that of the the ADS/CFT correspondence 176j which has been studied in the context of graphene 
in [l77l | . The search for charge fractionalization both for purely theoretical reasons as well as for potential appli cation 
in the field of quantum computing has also used graphene as a simple example to show the main feature s Il78ll . The 
same can be said on the search for Majorana fermions, another exploding subject in condensed matter |179l | with a 
clear relation with its high energy counterpart. 
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